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We use one-loop chiral perturbation theory to compare potential lattice computations of the — > 
tt^tt" decay amplitude at mK = 2m^ with the experimental value. We find that the combined 
one-loop efi'ect due to this unphysical pion to kaon mass ratio and typical finite volume effects is 
still of order minus 20-30%, and appears to dominate the effects from quenching. 



PACS numbers: f3.25.Es, f 2.38.Gc, 12.39. Fe 
I. INTRODUCTION 

Recently, we have used chiral perturbation theory 
(ChPT) to one loop in order to investigate three sys- 
tematic effects which affect lattice computations of the 
weak matrix element for K'^ — > tt+tt'^ decay: quench- 
ing, finite-volume effects, and the use of unphysical val- 
ues of the quark masses and pion external momenta |^] 
(to which we will refer as I). (For an extension to the 
partially quenched case, see ref. [||; for a recent lattice 
computation, see ||]; for other references, see /.) We con- 
sidered the case of a lattice computation (extrapolated to 
the continuum limit) with three degenerate light quarks, 
and final-state pions at rest. This is unphysical, because 
SU(3) flavor is broken in the real world, and because this 
choice of external momenta does not conserve energy. 

Here, we extend our results to the case mx = 2m7r, 
considering both the quenched and unquenched theories. 
This choice of masses conserves energy when all external 
mesons are at rest, and, since in reality mx ~ 3.6to^, 
it is closer to the real-world meson masses. One might 
therefore expect the systematic errors for this choice of 
masses to be smaller than in the energy-nonconserving 
case with degenerate masses. 

The choice mx = Smjr was also advertized in ref. [Q, 
where it is used as one of the ingredients in an improve- 
ment program for lattice computations of nonleptonic 
kaon decays with Wilson fermions. 

Here, we restrict ourselves to a summary of our results, 
and a brief discussion of differences with /; for other de- 
tails we refer to the extensive explanations contained in 
/. We also discuss similar systematic effects for Bk- The 
notation is the same as that of /. 



II. LATTICE METHOD 

On the lattice, usually operators with zero spatial mo- 
mentum, 0{t) — J2x^i^^^)' ^^'^ used. For — tt+tt*' 
decay on the lattice, one computes the time-correlation 
function 

C{t2,h)= {0\TT+it2)7T°it2)O4{h)K~{0)\0) (H.l) 

(0|^+(0)7r°(0)|^+^°)(j^+|i^-(0)|0) 
' (7r+7rO|7r+7rO)(i4:+|if+) ^ 



(7r+7r°|O4(0)|i^+)e-^^''(*^-*i)e-™^*i , (H.2) 

where O4 is the AS — 1, A/ — 3/2 part of the weak effec- 
tive hamiltonian Q . Eq. ( pl.2| ) shows the dominant term 
for large time separations, from which the desired K^- 
decay matrix element (7r+7r°|04|i^"'") can be extracted. 
This matrix element represents the decay process, con- 
fined to a finite lattice volume (we will restrict ourselves 
to a cubic volume with linear dimension L and periodic 
boundary conditions), in which a kaon, at rest, decays 
into a state consisting of two pions at rest with the low- 
est energy £'27r (s). The case for which niK — 2m.„ is 
of special interest, since, in that case, the matrix el- 
ement corresponds to an energy-conserving process (in 
the infinite- volume limit), in contrast to, for instance, 
the mass-degenerate case {rriK = ?n-n-) in which energy is 
injected through the weak operator O4. 

The aim of this paper is to compare real- world physical 
quantities with those obtained from a hypothetical lattice 
computation in full or quenched QCD with — ^ttIt,, 
using ChPT. Here, we will refer to the choice rriK = 2to^, 
as "lattice masses." Also, we will only consider the case 
with unbroken isospin, m„ — 7^ rris- 

III. ANALYTIC RESULTS FOR tuk = 2m^ FROM 
ONE-LOOP CHPT 

To 0(p4) in (quenched) ChPT, the one-loop diagrams 
(a) to (d) of Fig. 1 in /, along with the relevant wave- 
function renormalizations, as well as the tree-level con- 
tributions of 0{p^) weak operators Q, have to be evalu- 
ated in order to obtain the -decay matrix element. In 
/, where the mass-degenerate case was considered, spe- 
cial care was taken to accommodate the kinematic situ- 
ation in which energy is not conserved. In contrast, for 
rriK — 2rn^ and in infinite volume, standard Feynman 
diagram techniques can be used. In a finite volume, it 
was shown in / that (power-like) finite-volume correc- 
tions come exclusively from diagram (b) of Fig. 1, in 
which the final-state pions from the weak decay of the 
kaon rescatter. We have derived the finite-volume cor- 
rections for rriK 7^ TOtt, and, in particular, for — 2m.^ 
in the same way. 
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The weak decay operator O4 does not couple directly to 
the singlet meson rjo. However, for to„ = ^ m^, mix- 
ing between rjs and rjo occurs. As a result, the 77 two-point 
function inherents the "double pole" of the quenched the- 
ory 1^ , introducing dependence on the parameters d and 
a through 77 loops {S = TOq/ (247r^/^), with mo the singlet 
part of the 77' mass), which is not present in the mass- 
degenerate case. It turns out that these contributions are 
finite (independent of the cutoff), and, therefore, unam- 
biguously predicted in quenched ChPT. 

The one-loop result for the unquenched case follows di- 
rectly from Eqs. (43,44) (re-expressing the in terms 
of the bare decay constant / using Eq. (7) of/) by substi- 
tuting mx = 2m,r- In the quenched case (which in / was 
only considered for mx = titt), we obtain for ttik — 2TO7r, 
ignoring the contributions from 0{p*) operators. 



{TT+n'^\04\K+)'^ 



9X0^27 n 

— — = — ml^ X 
-2 log ^+CK + iG(m^L) 



(III.l) 
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with 



CK 



Note that the function G{x) is different from the one that 
appears in the mass-degenerate case, c/. Eq. (84) in /. 

For rriK 7^ Wtt, 0{p^) weak operators lead to contribu- 
tions to the one- loop correction factor of the form {Am'^ + 
BrriT^mK + Cm\)/{^TTfqf' . In the mass-degenerate 
limit, this leaves a dependence on the linear combina- 
tion A + B + C] whereas for ttik = Imj^ the linear 
combination is A + 2B + 4C. In Eq. ( 111.1 ), we have 
set A + 2B + AC — 0. Needless to say, for a comparison 
of complete 0{p*) results between different masses, infor- 
mation will be needed about the values of the coefhcients 
A, B and C. 



IV. COMPARISONS 

In this section, we will compare a hypothetical lattice 
computation, in a finite volume and with niK = 2.7X1^^^ 
of the TT+vr" matrix element and Bk with the 

real world. We will denote the lattice meson masses 
by niKMtt and m.,r,iatt = \mKMtt, and the real- world 
masses by and rriK, with — 136 MeV and 
TOj<- = 496 MeV. We wiU use f^^iatt = A = 132 MeV 
(/^ will only appear in one- loop corrections). The mass 
of the 77 is always determined from the tree-level relation 



to; 



K 



iTO^ 



Little is known about the values of 0(p*)-operator co- 
efficients. Therefore, as in /, we will ignore contributions 



— log 4 , 
12 ^ 12 



^ ^ emciems. increiore, as m j , we win ignore comriouiions 

— log? VSarctan — , (HI. 2) from 0{p'^) operators, and use the values 770 MeV and 



Q = -^-21og4+i^log7, 
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14 431 
c. = y + 21og4-— log7. 
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and 



1 GeV for the cutoffs A and A^ of the full and quenched 
theories respectively, and take the spread as an indication 
of the uncertainty from the lack of information about the 
0{p'^) constants. In subsection IV A, we will consider the 
case of a fuU-QCD lattice computation, and in subsection 
IVB the quenched case. 
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17.827 IOtt^ 
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The super/subscript q denotes quenched quantities: fq is 
the bare decay constant of the quenched theory, etc. The 
first factor of Eq. (III.l) is the tree- level result; the factor 
in parentheses gives the one-loop correction factor. 

Be fore substituting rrix = 2to^ in order to obtain 
Eq. (ni.l), one finds nonanalytic contributions coming 



A. Full QCD at rriK = 2m^ 

In the full theory, the 7r+7r° matrix element for 

the real world (subscript phys) and on the lattice (sub- 
script latt) can be related to one loop using Eqs. (43,44,7) 
of / (which was derived for arbitrary rriK and to^). For 
the real world, one evaluates the result for {■K^'K'^\Oi\K'^) 
at the physical values of rriK and to^t; whereas on the 



from poles at 



and 2to 



K 



mass of a pure ss meson). These nonanalytic contribu- 
tions have the values ck, cs and Cq after substituting 
rriK = 2m7r. The reason that we do not display the more 
general result here is that it is different depending on 
whether energy is conserved, or all external mesons are 
at rest. Only at rriK = 2m.„ do both more general ex- 
pressions agree with each other. 

The finite-volume correction (term proportional to 
G{inT^L) in Eq. (III.l)) has the same form in the full and 



(which is the lattice, one chooses m„ 



Jatt 

m- = 5to^ ;^tt/4. We obtain 



h'lTT'K latt and therefore 



{TT+7T°\0i\K 



/ phys 



X 



4 TO 



K 



^K.latt 



latt ' 
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where the superscript / denotes the full theory, and 
where 



quenched theories, since it comes only from the pion-pion 
rescattering diagram, which has no internal quark loops. 



X 



1 + a, 



phys 



+ Ulatt + 4(4^/Cy^ G{mK.lattL/2) 



(IV.2) 
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is the one-loop correction to the tree-level "conversion 
factor" 4(to|- — m^)/(3m|. ^^jj {U denotes the one-loop 
correction term (c/. Eq (87) of /) in infinite volume, 
i.e. with G = 0). We will restrict ourselves to the case 
where the lattice kaon mass is the same as the physical 
kaon mass, i.e. mK,iatt = Tnx, which will presumably be 
accessible in future lattice computations. It is straight- 
forward to consider other examples. 

The tree- level conversion factor is unambiguous 
and, for our example, equal to 1.23. We will therefore 
concentrate on the one-loop factor X. From Eq. (87) of /, 
Uphys = 0.0888 and -0.0146 for A = 1 GeV and 770 MeV 
respectively (we will ignore the imaginary part of the ma- 
trix element, since it does not contribute to the magni- 
tude of the amplitude to order p^). For ■mK,iatt = niK, 
we obtain Uiatt = 0.328 and 0.147 for A = 1 GeV and 
770 MeV respectively. 

On the lattice, when L is such that utikL = 6 or 8, the 
relative one-loop contribution from the finite volume cor- 
rection 



G{mKL/2) is 0.215 or 0.134, respectively. 



Since, as explained before, we have omitted different lin- 
ear combinations of O(p^) coefficients in the real-world 
and lattice cases respectively, we will vary the values of 
the cutoff in Uphys and Uiatt independently. 

Wc list in Table | the values of X for four combinations 
of the values 1 GeV or 770 MeV for the cutoff, and for 
volumes such that itlkL = 6, ttikL = 8 and itlkL = oo. 
We take the spread of the factor X due to changes in the 
values of the cutoff as a systematic error, which, from 
Table |, is around 15 - 20%. 

The one-loop expression for Bk in the full the- 
ory is given by Eq. (36) of /. For physical masses. 
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f, phys 
K 



-- 1 GeV or 770 MeV, 



/Bi = 1.72 or 1.42 for A = 
respectively; whereas for mK,iatt 
^Uatt^gf ^ ^ .,3 ^^^^^ j^^, ^j^g gg^^g values of the 

cutoff. If we compare at the same value of the cutoff, we 
see that the real-world and lattice values differ by about 
1%, in contrast to the if^-decay matrix element. How- 
ever, if we compare at different values of the cutoff, again 
as an estimate of the error introduced by ignoring O(p^) 
coefficients, we see that the real-world and lattice values 
may differ by as much as 20%. 

Following /, we can examine the ratio 

\fK {7r+7r°\0,\K+)/ {k"\0'\K°)Y 

J latt 



(IV.3) 



rriKL 


^(1) 


„{0.77) 
^(0.77) 


^(0.77) 


y(0.77) 


6 


0.71 


0.72 


0.80 


0.64 


8 


0.74 


0.77 


0.85 


0.67 


00 


0.82 


0.86 


0.95 


0.74 



TABLE I. The factor X for different values of itikL 
and different combinations of values of the cutoff A (su- 
per/subscripts (1) and (0.77) denote cutoff used in the nu- 
merator/denominator of Eq. ([V.2)). 



which is independent of the 0(p^)-operator coefficient 
a27. The tree-level value for the above ratio is 9z/8\/2- 
At one loop, we find, for mx.iatt = '^rriTr^iatt = mx, cor- 
rections of -53%, -61% and -74% for A = 1 GeV and 
TO^L = 6, 8 and 00, respectively (for A = 770 MeV, the 
corrections are —29%, —37% and —50%). 



B. Quenched QCD at tuk ~ 2m,r 

We now compare real-world quantities and quantities 
as would be obtained from a quenched lattice computa- 
tion with rriK = 2mT^. We get, from Eq. (III.l), 



{TT+AOi\K+)'i 



unphys 



9ia27 



(IV.4) 
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'^latt ~ 



mK,latt 
WW 



-21og- 



'■K,latt 
"A2~ 



0.0697(5 + 0.0603a 



0.6820 



''Kdatt 
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In analogy with Eq. ( IV.l ), we relate the real- world 
decay matrix element to that of the quenched lattice 
theory, obtaining 



{n+n"\0,\K+)i 
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where 



Y = 



1 + U, 



phys 



^ + uLt 



4(4^/°"^ G{mK.lattL/2) 
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is again the one-loop correction to the tree-level conver- 
sion factor. 

Again, we will restrict ourselves to the case mx.iatt — 
'^n^TT^iatt = TTiK, and take fq — fj^. First, we estimate 
the imp ortan ce of the 5- and a-dependent terms in Uf^^^^ 
cf. Eq. ( [V.5 ). The S- and a- independent part of U^^^^ is 
equal to 0.69 and 0.59 for A, = 1 GeV and A, = 770 MeV 
respectively. The value for S is estimated to be less than 
around 0.2 |jl^. a is only poorly known, but is unlikely 
to be larger than one in magnitude [|l^. Therefore, the 
total contribution of the S and a terms is less than about 
ten percent of the S- and a-independent contribution. 
We will omit the S and a terms in Table || below. Table 
U gives the values of Y for the four combinations of the 
values 1 GeV or 770 MeV for the cutoffs, and for volumes 
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such that mifl/ = 6, itikL = 8 and itikL = oo. These 
values of Y deviate substantially from the tree-level value 

Y — 1, like in the examples with niK = TOtt considered 
in /. We see that the spread in the values of the factor 

Y is 15% or less. 

For the quenched one- loop expression for Bk , we refer 
to Eq. (37) of/. For mK.iatt = '2mTrjatt = rriK, 



B 



q.latt 
K 



Bi 



1-0.14(5 



(IV.8) 



0.691 
0.358- 



0.050Q, 
0.085a, 



Kq = 1 GeV 
A, = 770 MeV 



The 5 and a terms are again relatively small. If we com- 
pare with the values of Bff'^^^/B-f (see previous subsec- 
tion), we see that again real- world and quenched lattice 
values differ only by a few percent if we compare at the 
same value of the cutoffs. At different values of the cut- 
offs, they may again differ by as much as 20%. 

In the quenched case, the ratio TZ defined in Eq. ( IV.3| ) 
is, for mK,iatt = '^m^Mtt = ™K and fq = 



7^o 



^ (7r+7r"\04\K^ 

JK =0 

(K \0'\K^>) 



latt 



9i 



(IV.9) 



1 + 0.056(5 + 0.0224G(mifL/2) 



-0.380 
-0.140- 



0.058ck, 
0.093ck, 



A, 
A„ 



1 GeV 
770 MeV 



We see that the quenched values of TZ are closer to the 
tree-level value than the unquenched values, for the same 
choice of meson masses and volumes. 



masses, and quenching (Y) for this matrix element. (The 
tree-level conversion factor corrects only for unphysical 
masses.) For more discussion on the reliability of such 
estimates, and other systematic errors, see /. 

In Tables | and we give numerical examples, illus- 
trating these results for a lattice kaon mass equal to the 
experimental value, rriK = 496 MeV. The four different 
values on each line represent four combinations of cut- 
offs we chose in evaluating the ratios X and Y, and we 
take the spread as an indication of the uncertainties in- 
troduced by the lack of information on O(p^) constants. 
We should emphasize that ChPT does not give us any in- 
formation about the ratio {0127 fq) / [alj p) in Eq. ( [V.6| ). 

From the Tables, we conclude that even in the "more 
physical" case mx = 'imj^, the one-loop systematic effect 
may still be as large as minus 20-30%. Comparing Ta- 
bles Q and we see that quenching seems to have only 
a minor effect. The large deviations from one of the ra- 
tios X and Y are caused by the sensitivity of the matrix 
element to the ratio mT^/rriK and finite- volume effects. 

This situation is different from that of Bk- Here, if 
we compare values of Bk between the real world and 
the lattice at the same value of the cutoff, the difference 
is only a few percent, both quenched and unquenched. 
However, comparing at different values, again in order to 
get an idea of the effect of the O(p^) constants, indicates 
that also here systematic effects can be as large as 20%. 

We thank Tanmoy Bhattacharya and Rajan Gupta for 
discussions. This work is supported in part by the US 
Department of Energy (MG by the Outstanding Junior 
Investigator program) . 



V. CONCLUSION 

Ou r main r esults are the ratios X and Y in 
Eqs. ( |IV.2| , |IV.7D , which give the one-loop ChPT correc- 
tion to the conversion factor between a lattice computa- 
tion of the 7^+ — !■ vr -|- TT*' amplitude with ttik — 2m^ and 
the experimental value, for unquenched and quenched 
QCD, respectively. These ratios give an estimate of the 
systematic effect due to finite volume, unphysical quark 







v(0-77) 
-^(0.77) 


y(l) 
(0.77) 


^(1) 


itikL = 6 


0.71 


0.69 


0.76 


0.65 


rnxL — 8 


0.75 


0.73 


0.81 


0.68 


00 volume 


0.83 


0.81 


0.89 


0.75 



TABLE II. The factor Y for different values of rriKL and 
different combinations of values of the cutoffs A and A^ (su- 
per/subscripts (1) and (0.77) denote cutoff used in the nu- 
merator/denominator of Eq. ([V.7)). 
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